clear
use "H:\MOUM2012\MOUM2012.dta"

global controlvars  lncprt  majpow cntgd dist  numstate  peaceyrs _spline1 _spline2 _spline3
**********************************************************************************
******************************AIC AND BIC TESTS***********************************
**********************************************************************************
logit crisis1 dml dmdist  $controlvars, cl(dyadid) nolog
estat ic
estat ic, n(13538)
fitstat
logit crisis1 bdm dmdist  $controlvars, cl(dyadid) nolog
estat ic
estat ic, n(13538)
fitstat
logit crisis1 dm10 dmdist  $controlvars, cl(dyadid) asis nolog
estat ic
estat ic, n(13538)
fitstat
logit crisis1 dmlsq dmdist  $controlvars, cl(dyadid) nolog
estat ic
estat ic, n(13538)
fitstat

**********************************************************************************
******************TABLE 1: CIE, DEMOCRACY AND ICB ONSET, 1960-2001****************
**********************************************************************************
quietly {
eststo clear
*Table 1, Model 1
eststo: logit 		crisis1 	   		  dml   dmdist 			  $controlvars if CIEl!=., cl(dyadid) nolog
*Table 1, Model 2
eststo: logit 		crisis1 		CIEl  dml   dmdist 			  $controlvars, cl(dyadid) nolog
*Table 1, Model 3
eststo: logit 		crisis1 		CIEl  bdm   dmdist			  $controlvars, cl(dyadid) nolog
*Table 1, Model 4
eststo: logit 		crisis1 		CIEl  dm10  dmdist			  $controlvars, cl(dyadid) nolog
*Table 1, Model 5
eststo: logit 		crisis1         CIEl  dmlsq	dmdist			  $controlvars, cl(dyadid) nolog
*Table 1, Model 6
eststo: logit 		crisis1 		CIEl  	                      $controlvars, cl(dyadid) nolog
}
esttab, b(2) se(2) replace label star(* 0.10 ** 0.05 *** 0.01)  order(CIEl dml bdm dmlsq dmdist $controlvars) scalars("ll Log lik.") pr2 varwidth(25) modelwidth(8)

**********************************************************************************
******************TABLE 2: CIE, DPL, CAPOPENL, PUBH AND ICB ONSET, 1960-2001******
**********************************************************************************
quietly {
eststo clear
*Table 2, Model 1
eststo: logit crisis1 CIEl 		 edvl 					$controlvars, cl(dyadid) nolog
*Table 2, Model 2
eststo: logit crisis1 CIEl 		 edvl	dpl			    $controlvars, cl(dyadid) nolog
*Table 2, Model 3
eststo: logit crisis1 CIEl 		 edvl 	capopenl_ipol2	$controlvars, cl(dyadid) nolog
*Table 2, Model 4
eststo: logit crisis1 CIEl 		 edvl 	pubh			$controlvars, cl(dyadid) nolog
}
esttab, wide b(2) se(2) replace label star(* 0.10 ** 0.05 *** 0.01)  order(CIEl edvl dpl capopenl_ipol2 pubh $controlvars) scalars("ll Log lik.") pr2 varwidth(25) modelwidth(8)

**********************************************************************************
******************TABLE A1: RESULTS WITHOUT CIEL IN TABLE 1***********************
**********************************************************************************
quietly {
eststo clear
*Table A1, Model 1
eststo: logit 		crisis1	    bdm	 	     dmdist	  $controlvars if CIEl!=., cl(dyadid) nolog
*Table A1, Model 2
eststo: logit       crisis1     dmlsq		 dmdist	  $controlvars if CIEl!=., cl(dyadid) nolog
}
esttab, wide b(2) se(2) replace label star(* 0.10 ** 0.05 *** 0.01)  order(bdm dmlsq dmdist $controlvars) scalars("ll Log lik.") pr2 varwidth(25) modelwidth(8)

**********************************************************************************
******************TABLE A2: RESULTS WITHOUT CIEL IN TABLE 2***********************
**********************************************************************************
quietly {
eststo clear
*Table A2, Model 1
eststo: logit crisis1  		 edvl 					$controlvars if CIEl!=., cl(dyadid) nolog
*Table A2, Model 2
eststo: logit crisis1  		 edvl	dpl			    $controlvars if CIEl!=., cl(dyadid) nolog
*Table A2, Model 2
eststo: logit crisis1  		 edvl 	capopenl_ipol2	$controlvars if CIEl!=., cl(dyadid) nolog
*Table A2, Model 2
eststo: logit crisis1  		 edvl 	pubh			$controlvars if CIEl!=., cl(dyadid) nolog
}

esttab, wide b(2) se(2) replace label star(* 0.10 ** 0.05 *** 0.01)  order(edvl dpl capopenl_ipol2 pubh $controlvars) scalars("ll Log lik.") pr2 varwidth(25) modelwidth(8)



